Sub-eV Hubble scale inflation within gauge mediated supersymmetry breaking 
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Minimal Supersymmetric Standard Model with gauge mediated supersymmetry breaking has all 
the necessary ingredients for a successful sub-eV Hubble scale inflation Hint ~ 10""^ — 10~^ eV. The 
model generates the right amplitude for scalar density perturbations and a spectral tilt within the 
range, 0.90 < < 1. The reheat temperature is Tr <^ 10 TeV, which strongly prefers electroweak 
baryogenesis and creates the right abundance of gravitinos with a mass 7713/2 ^ 100 keV to be the 
dark matter. 



Embedding a very low scale inflation within a particle 
physics model is a challenging problem. It is not only 
difficult to obtain sufficient number of e- foldings, right 
amplitude for the scalar density perturbations, the right 
tilt in the power spectrum, but also generating baryon 
asymmetry and dark matter simultaneously. 

Recently there has been a real progress in our under- 
standing of embedding inflation within particle physics, 
particularly within the Minimal Supersymmetric Stan- 
dard Model (MSSM), where the inflaton belongs to the 
MSSM instead of being an ad-hoc gauge singlet. Fore- 
most the models 0, not only predict the right ampli- 
tude of the scalar density perturbations and the tilted 
spectrum, see also Q , but are also testable at LHC ^ . 

In this paper we provide a simple example of a sub- 
eV Hubble scale inflation. Hint ~ 10""^ — 10~^ eV, em- 
bedded within MSSM. This is realizable provided super- 
symmetry (SUSY) breaking is communicated via gauge 
mediation, i.e. gauge mediated supersymmetry breaking 
(GMSB) H, in contrast to Refs. 0,0 where we assumed 
gravity mediation. We shall predict the reheat tem- 
perature around TeV, which strongly favors electroweak 
baryogenesis within MSSM and gravitino as a dark mat- 
ter candidate. Thus all the ingredients for a successful 
cosmology are naturally contained within MSSM. 

Let us first highlight relevant points of the model: 

• The model is based on inflation in the vicinity of 
a saddle point, see 0,0, which is predictive, ra- 
diatively stable, free from supergravity and Trans- 
Planckian corrections. The only scales which enters 
in the potential is the weak scale ^ TeV, and the 
Grand Unified scale, Mgut lO^^'^ GeV. 



^ The inflaton candidates are, LLe, udd [Jj and NHuL Q flat 
directions. Here L denotes left-handed sleptons, u, d, e denote 
right-handed squarks and sleptons, respectively, TV denotes right- 
handed sneutrinos and Hu is the Higgs which gives masses to the 
up type quarks. The LLe, udd mass should be > 340 GeV 0, 
which is within the reach of LHC. For other attempts of infla- 
tion with flat directions, see |^|^. For a review on MSSM flat 
directions, see 0. 



• The reheat temperature of the model is just around 
TeV, which points towards thermal electroweak 
baryogenesis y. If the reheat temperature is 
slightly lower than this, we have the right condi- 
tions for a cold electroweak baryogenesis 10]. 

• Within GMSB the gravitino is the hghtest SUSY 
particle (LSP). Therefore a reheat temperature of 
©(TeV) wiU lead to a sufficient relic abundance 
for the (stable) gravitino as a dark matter candi- 
date [U. 

Let us now consider an MSSM flat direction, (j), lifted 
by a non-renormalizable {n > 3) superpotential term (for 
a detailed dynamics on multiple flat direction, see [ia|): 

W = {XJn){<i>yM-^^), (1) 

where $ is the superfield which contains the flat direction 
(j). Within MSSM all the flat directions are lifted by n < 9 
operators The cut-off scale is Mgut, therefore 

the above superpotential is a reflection of integrating out 
the physics above the GUT scale, and we assume the non- 
renormalizable coupling to be A„ ~ 0{1). Note that the 
new physics does not necessarily have to be tied to the 
GUT physics. For example a gauged U{1)b-l may ap- 
pear at an intermediate scale Mb-l «C Mqut- This will 
amount to the same superpotential as in Eq. Q parame- 
terized by Mgut, so long as A„ ^ 1. This is conceivable 
as A„ is a product of the Yukawa couplings associated 
with the new interaction terms beyond the MSSM. 

In GMSB the two-loop correction to the flat direction 
potential results in a logarithmic term above the messen- 
ger scale, i.e. > Mg |0. Together with the ^-term 
this leads to the scalar potential 

V - Ml In [Jf2)+Aco.{r^e+e^)-^+Xl^^ , 

V^"S/ "^^^GUT ^^GUT 

(2) 

where Mp ~ [msuSY x Mg)^/^ and msusY ~ 1 TeV 
is the soft SUSY breaking mass at the weak scale. For 
(f) > A'fp/ 171^/2, usually the gravity mediated contribu- 
tion, 7713/2'/'^, dominates the potential where 7713/2 is the 
gravitino mass. Here we will concentrate on the VEVs 
M, < < Mj./ 1713/2. 
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In Eq. ||2Jl, (f> and 9 denote the radial and the angular 
coordinates of the complex scalar field $ = exp[i9] re- 
spectively, while 9a is the phase of A-teim (thus A is a 
positive quantity with a dimension of mass). Note that 
the first and third terms in Eq. ||2J) are positive definite, 
while the A-tevm leads to a negative contribution along 
the directions where cos{n9 + 9a) < 0. The cosniological 
importance of an A-teim can be found in (ij, |l^ . 

Although individual terms are unable to support a sub- 
Planckian VEV infiation, but as shown in Refs. Q, H, H, 
01, a successful inflation can be obtained near the saddle 
point, which we find by solving, V'{(j)o) — F"(0o) = 
(where derivative is w.r.t 0). 



Af "-3 Af 2 
GUT ^ 



A = 



4(n- 1)2A„ ,„_2 



i/(«-i) 



(3) 
(4) 



In the vicinity of the saddle point, we obtain the total 
energy density and the third derivative of the potential 
to be: 



In 



Ml 



3n- 2 



n(n — 1) 



V"'{cbo) = 4n(n-l)M|0o = 



(5) 
(6) 



There are couple of interesting points, first of all note 
that the scale of inflation is extremely low in our case, 
barring some small coefficients of order one, the Hubble 
scale during infiation is given by: 



Hi, 



10-'^ - 10^^ eV. 



(7) 



for Mp '--^ 1 — 10 TeV. For such a low scale infiation usu- 
ally it is extremely hard to obtain the right phenomenol- 
ogy. But there are obvious advantages of having a low 
scale infiation. Alp ^ -ffinf- The supergravity corrections 
and the Trans-Planckian corrections are all negligible 4] , 
therefore the model predictions are trustworthy. 

Perturbations which are relevant for the COBE nor- 
malization are generated a number A/cobe e-foldings be- 
fore the end of inflation. The value of A/cobe depends on 
thermal history of the universe and the total energy den- 
sity stored in the inflaton, which in our case is bounded 
by, Vo < 10^^ (GeV)'*. The required number of e-foldings 
yields in our case, A/cobe ~ 40 0|, provided the uni- 
verse thermalizes within one Hubble time. Although 
within SUSY thermalization time scale is typically very 
long ,19j . however, in this particular case it is possible to 
obtain a rapid thermalization. 

Near the vicinity of the saddle point, (j)o, the poten- 
tial is extremely flat and one enters a regime of self- 
reproduction 20] . The sclf-rcproduction regime lasts as 
long as the quantum diffusion is stronger than the clas- 
sical drag; i/inf/27r > (p/Hini, for 0s < (t> ^ 4>o, where 
(^0 — f^s — Mp {(po/Mp)'^^^ . From then on, the evolution 



is governed by the classical slow roll. Inflation ends when 
\r]\ 1, which happens at ~ 0e, where 



PO = 



2n{n- l)Mp 



(8) 



Assuming that the classical motion is due to the third 
derivative of the potential, V'{(f>) ~ {1/2)V"' {(f>o){(j) - 
(/)o)^, the total number of e-foldings during the slow roll 
period is found to be: 

J^tot — / ■■ — -; — 7 

This simplifles to 



2Vocj)l 



^08 ( 1 

l)Mj,M^ \(l)o- 



(9) 
(10) 



Let us now consider the adiabatic density perturbations. 
Despite Hi^i ^ 1 eV, the fiat direction can generate ad- 
equate density perturbations as required to explain to 
match the observations. Recall that infiation is driven 
by V" ^ 0, we obtain 



(l/5^)(iJ?f/^) ~ M|MpAA20o' 



10" 



(11) 



Note that for Mp ~ 10 T eV, and A/cobe ~ 40, we match 
the current observations [2l|, when (j)o '-^ 10^^ GeV. The 
validity of Eq. Q for such a large VEV requires that 
Mj, > (10" GeV) X m3/2. For Mp - 10 TeV this yields 
the bound on the gravitino mass, < 1 MeV, which 

is compatible with the dark matter constraints as we will 
see. 

We can naturally satisfy Eq. (|ll|l provided, n — 6. The 
non-renormalizable operator, n = 6, points towards two 
MSSM flat directions out of many. 



LLe and 



udd . 



(12) 



As we discussed before in .1] , these are the only directions 
which are suitable for inflation as they give rise to a non- 
vanishing A-term. Note that the inflatons are now the 
gauge invariant objects. The total number of e-foldings, 
during the slow roll infiation, after using Eq. (|10|) yields. 



Kot - 10^ 



(13) 



While the spectral tilt and the running of the power spec- 
trum are determined by A/cobe ~ 40 <C Atot ■ 



= 1 + 2?] — 6e ~ 1 



0.90, 



COBE 



(14) 



drig 
dink 



16677 - 246^ - 2C 



-'^cobe 



-10"^, 



(15) 



where = Mj,V'V"' /V^. Note that the spectral tilt is 
slightly away from the 2a result of the current WMAP 3 



3 



years data, on the other hand running of the spectrum is 
well inside the current bounds [2ll |. 

At first instance one would discard the model just from 
the slight mismatch in the spectral tilt from the current 
observations. However note that our analysis strictly as- 
sumes that the slow roll inflation is driven by V"'{(j)o). 
This is particularly correct if y (0o) — and V"{4)a) = 0. 
Let us then study the case when V'{(f>o) ^ 0, as discussed 
in 0. 

The latter case can be studied by parameterizing a 
small deviation from the exact saddle point condition by 
solving near the point of inflection, where we wish to solve 
V"{(t)o) — and we get up-to 1st order in the deviation, 
(5 < 1, 



ijo = <po \ ^ —S 



n{n - 2) 



(16) 



with A and <j>Q are the saddle point solutions. Then the 
1st derivative is given by 



(17) 



Therefore the slope of the potential is determined by, 
y (</>) :^ V'iM + (l/2)l^"'((/.o)(</> - M'- 

Note that both the terms on the right-hand side are 
positive. The fact that V'{(j)o) ^ can lead to an in- 
teresting changes from the saddle point behavior, for in- 
stance the total number of e-foldings is now given by 



Mot = 



d(j) 



(18) 



First of all note that by including V' , we are slightly 
away from the saddle point and rather close to the point 
of inflection. This affects the total number of e-foldings 
during the slow roll. It is now much less than that of 
Mot, i.e. Mot < 10^ see Eq. (O. 

When both the terms in the denominator of the inte- 
grand contributes equally then there exists an interesting 
window. 



- < 6 < - . 
8 - - 2 



(19) 



where 



n{n — ly 



In 



Ml 



3n- 2 

n(n — 1) 



^^pMaoBE 



(20) 
0, 



The lower limit in Eq. (|19|l is saturated when T^'((/)o) 
while the upper limit is saturated when Mtot — Mcobe — 
40. It is also easy to check that there will be no self- 
reproduction regime for the field values determined by 
5. 

It is a straightforward but a tedious exercise to demon- 
strate that when the upper limit of Eq. (|19|l is satu- 
rated the spectral tilt becomes ~ 1, when the lower 
limit is satisfied we recover the previous result with 



Us = 0.90. This value, rig — s- 1, can be easily understood 
as (/'coBE 4>Q (where (/icobe corresponds to the VeV 
where the end of inflation corresponds to Mcobe ^ 40), 
in which case, 77 ^ 0. Therefore the spectral tilt becomes 
nearly scale invariant. We therefore find a range 



0.90 < < 1 



(21) 



whose width is within the 2a error of the central 
limit [231 ■ Similarly the running of the spectral tilt gets 
modified too but remains within the observable limit ^, 
while the amplitude of the power spectrum is least af- 
fected Jj]. 

Let us now discuss the issue of reheating and thermal- 
ization. Important point is to realize that the inflaton be- 
longs to the MSSM, i.e. LLe and udd, both carry MSSM 
charges and both have gauge couplings to gauge bosons 
and gauginos. After inflation the condensate starts os- 
cillating. The effective frequency of the inflaton oscil- 
lations in the Logarithmic potential, Eq. 0, is of the 
order of Mp/(pQ, while the expansion rate is given by 
Hint ~ Mp/Mp. This means that within one Hubble 
time the inflaton oscillates nearly Mp/(j>o ^ 10^ times. 
The motion of the inflaton is strictly one dimensional 
from the very beginning. During inflation, the imaginary 
direction is very heavy and settles down in the minimum 
of the potential. 

An efficient bout of particle creation occurs when the 
inflaton crosses the origin, which happens twice in ev- 
ery oscillation. The reason is that the fields which are 
coupled to the inflaton are massless near the point of en- 
hanced symmetry. Mainly electroweak gauge fields and 
gauginos are then created as they have the largest cou- 
pling to the flat direction. The production takes place in 
a short interval. Once the inflaton has passed by the ori- 
gin, the gauge bosons/gauginos become heavy by virtue 
of VeV dependent masses and they eventually decay into 
particles sparticles, which creates the relativistic ther- 
mal bath. This is so-called instant preheating mecha- 
nism j22|. In a favorable condition, the flat direction 
VeV coupled very weakly to the flat direction inflaton 
could also enhance the perturbative decay rate of the in- 
flaton [231 . In any case there is no non-thermal gravitino 
production as the energy density stored in the infla- 
ton oscillations is too low. 

A full thermal equilibrium is reached when o) kinetic 
and b) chemical equilibrium are established |l9l |. The 
maximum temperature of the plasma is given by 



1/4 



Mf < lOTeV, 



(22) 



when the flat direction, either LLe or udd evaporates 
completely. This naturally happens at the weak scale. 



A similar exercise can be done for the running of the spectral 
tilt and the running lies between — IG/A/qq^j; < dus/dlnk < 



4 



There are two very important consequences which we 
summarize below. 

Hot or cold electroweak Baryogenesis: The model 
strongly favors electroweak baryogenesis within MSSM. 
Note that the reheat temperature is sufficient enough for 
a thermal electroweak baryogenesis . 

However, if the thermal electroweak baryogenesis is 
not triggered, then cold electroweak baryogenesis is still 
an option [l^. During the cold electroweak baryogene- 
sis, the large gauge field fluctuations give rise to a non- 
thermal sphaleron transition. In our case it is possible 
to excite the gauge fields of SU{2)l x U{\)y during in- 
stant preheating provided the inflaton is LLe. The LLe 
as an inflaton carries the same quantum number which 
has a B — L anomaly and large gauge field excitations 
can lead to non-thermal sphaleron transition to facilitate 
baryogenesis within MSSM. 

Gravitino dark matter. Within GMSB gravitinos are 
the LSP and if the i?-parity is conserved then they are 
an excellent candidate for the dark matter. There are 
various sources of gravitino production in the early uni- 
verse I25I I . However in our case the thermal produc- 
tion is the dominant one and mainly helicity ±1/2 grav- 
itinos are created. Gravitinos thus produced have the 
correct dark matter abundance for j3 US 

"^^/^ ~_ U ]( Y .23^ 

100 keV few VI TeV/ VI TeV/ ' ^ ' 

where Mg is the gluino mass. For TO3/2 ^ 100 keV, 



Eq. ^ is easily satisfied for Mg - 1 TeV and Tr < 10 
TeV. We remind that for O(keV) < 771.3/2 < lOolkeV 
gravitinos produced from the sferniion decays overdose 
the universe |15| |. 

Before concluding we should also highlight that the 
existence of a saddle point does not get spoiled through 
radiative corrections, see Ref. To summarize, we pro- 
vided a truly low scale inflation model with Hini ~ 10~^— 
10-1 eV, embedded within MSSM, provided GMSB is the 
correct paradigm. Although inflation occurs at such low 
scales, the model predictions match the current WMAP 
data and the reheat temperature of Tr <^ 10 TeV is suf- 
ficient enough to trigger either hot or cold electroweak 
baryogenesis. The model also produces sufficient abun- 
dance of gravitinos to be the dark matter candidate. 
Thus inflation within GMSB connects the physics of mi- 
crowave background radiation to a successful dark mat- 
ter and a baryogenesis scenario whose ingredients are 
testable at the LHC. 

We wish to thank Alex Kusenko and Misha Sha- 
poshnikov for helpful discussion. The research of RA 
was supported by Perimeter Institute for Theoretical 
Physics. Research at Perimeter Institute is supported 
in part by the Government of Canada through NSERC 
and by the provine of Ontario through MEDT. The re- 
search of AM is partly supported by the European Union 
through Marie Curie Research and Training Network 
"UNIVERSENET" (MRTN-CT-2006-035863) . 



[1] R. AUahverdi, K. Enqvist, J. Garcia-Bellido and 

A. Mazumdar, arXiv:hep-ph/0605035 
[2] R. AUahverdi, A. Kusenko and A. Mazumdar, 

arXiv:hep-ph/0608138 
[3] R. AUahverdi and A. Mazumdar, arXiv:hep-ph/0610069 
[4] R. AUahverdi, K. Enqvist, J. Garcia-Bellido, A. Jokinen 

and A. Mazumdar, arXiv:hep-ph/0610134 
[5] G. Lazarides and Q. Shafi, Phys. Lett. B 308, 17 (1993). 

S. Kasuya, T. Moroi and F. Takahashi, Phys. Lett. 

B 593, 33 (2004). R. Brandenberger, P. M. Ho and 

H. C. Kao, JCAP 0411, Oil (2004). A. Jokinen and 

A. Mazumdar, Phys. Lett. B 597, 222 (2004). 
[6] J. C. B. Sanchez, K. Dimopoulos and D. H. Lyth, 

arXiv:hep-ph/0608299 
[7] For reviews, see K. Enqvist and A. Mazumdar, Phys. 

Rept. 380, 99 (2003); M. Dine and A. Kusenko, Rev. 

Mod. Phys. 76, 1 (2004). 
[8] G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 

(1999). 

[9] V. A. Rubakov and M. E. Shaposhnikov, Usp. Fiz. Nauk 
166, 493 (1996) [Phys. Usp. 39, 461 (1996)]; 

[10] J. Garcia-Bellido, D. Y. Grigoriev, A. Kusenko and 
M. E. Shaposhnikov, Phys. Rev. D 60 (1999) 123504. 

[11] F. D. Steffen, arXiv:hep-ph/0605306 

[12] K. Enqvist, A. Jokinen and A. Mazumdar, JCAP 0401, 
008 (2004). 

[13] M. Dine, L. Randall and S. Thomas, Phys. Rev. Lett. 75, 



398 (1995). M. Dine, L. Randall and S. Thomas, Nucl. 

Phys. B 458, 291 (1996). 
[14] T. Gherghetta, C. Kolda and S. P. Martin, Nucl. Phys. 

B 468, 37 (1996). 
[15] A. de Gouvea, T. Moroi and H. Murayama, Phys. Rev. 

D 56, 1281 (1997). 
[16] M. Postma and A. Mazumdar, JCAP 0401, 005 (2004). 
[17] R. AUahverdi, K. Enqvist, A. Jokinen and A. Mazumdar, 

arXiv:hep-ph/0603255 
[18] C. P. Burgess, R. Easther, A. Mazumdar, D. F. Mota 

and T. Muhamaki, JHEP 0505, 067 (2005). 
[19] R. AUahverdi and A. Mazumdar, arXiv:hep-ph/0505050 

R. AUahverdi and A. Mazumdar, arXiv:hep-ph/0512227 

R. AUahverdi and A. Mazumdar, arXiv:hep-ph/0603244 

R. AUahverdi and A. Mazumdar, arXiv:hep-ph/0608296 
[20] A. D. Linde, PARTICLE PHYSICS AND INFLATION- 
ARY COSMOLOGY (Harwood Academic Publishers, 

Chur, Switzerland 1990). 
[21] D.N. Spergel, et.al., astro-ph/0603449 
[22] G. N. Felder, L. Kofman and A. D. Linde, Phys. Rev. D 

59, 123523 (1999). 
[23] R. AUahverdi, R. Brandenberger and A. Mazumdar, 

Phys. Rev. D 70, 083535 (2004) arXiv:hep-ph/0407230 '. 
[24] A. L. Maroto and A. Mazumdar, Phys. Rev. Lett. 84, 

1655 (2000). 

[25] R. AUahverdi, S. Hannestad, A. Jokinen, A. Mazumdar 
and S. Pascoh, |arXiv:hep-ph/0504102| 



M. Bolz, A. Brandenburg and W. Buchniiillcr, Nucl. 
Phys. B 606, 518 (2001). 



